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Abstract. We present a review of the standard paradigm for giant planet 
formation, the core accretion theory. After an overview of the basic concepts of 
this model, results of the original implementation are discussed. Then, recent 
improvements and extensions, like the inclusion of planetary migration and the 
resulting effects are discussed. It is shown that these improvement solve the 
"timescale problem" . Finally, it is shown that by means of generating synthetic 
populations of (extrasolar) planets, core accretion models are able to reproduce 
in a statistically significant way the actually observed planetary population. 



1. Introduction 

Our current understanding of planet formation is based on several centuries 
of observations of the planets of our own Solar System, 12 years of extrasolar 
planets detection, and several decades of observations of young stellar systems. 
These studies have let to the general concept that after the collapse of a dense gas 
cloud, a protostar surrounded by a protoplanetary disk was formed. In this disk, 
solids started to coagulate from fine dust and grew further by mutual collision to 
form planetesimals (provided the bottleneck by bodies roughly one meter in size 
can be overcome), then protoplanets, and ultimately the actual planets. Some 
of the protoplanets managed to accrete a massive gaseous envelope onto their 
core. This is the very rough outline of the core accretion model. 



2. Core Accretion Paradigm 

In the core accretion model, the formation of gas giant planets is therefore seen 
as a two step process: In the first one a solid core is formed which, provided 
this core reaches critical mass, gives rise to the second step in which runaway 
gas accretion occurs, leading to a quick build-up of a massive envelope. This 



basic scenario was first studied more than 30 years ago (jPerri &: Cameronlll974l ; 
Mizuno et aT~lll978l : iBodenheimer k Pollack! ll986T ) 



2.1. Baseline Formation Models 

The growth of the core (at a fixed semimajor axis in baseline models) occurs 
through collisional accretion of background planetesimals, which themselves are 
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Figure 1. Time evolution of the mass of accreted planetesimals (solid line), 
the mass of accreted gas (dotted line) and total mass of the planet (dashed 
line), for init ial conditions e quiva lent to the preferred model f or the formation 
of Jupiter in lPollack et all (|1996l ) (from I Alibert et al.ll2005ah . 



formed by collisional coagulation of small dust grains jWetherill fc Stewadfl989h 
or an instability in the dust layer (jjohansen et alj|2007l). The g o verni ng equation 
for the growth rate is given by the classical theory of Safronov ( 19691 ). However, 
in many early studies, dM core /dt was treated as a free parameter. The formation 
of the core occurs through the same mechanism as the one generally accepted 
for the formation of terrestrial planets. In this sense, core accretion is an unified 
theory for the formation of both terrestrial and giant planets. When the core 
has reached roughly the mass of the moon, it can hold an initially tenuous hy- 
drostatic atmosphere. Its structure is well described by the classical set of ID 
stellar structure equations except for the nuclear energy release term that has 
to be replaced by heating due to infalling planetesimals. To solve the structure 
equations, boundary values must be specified. In early models the outer bound- 
ary conditions, the background nebula values, were often taken to be constant, 
and opacities similar to the one o f the interstellar med ium were used. 

The calculations presented in Pollack et al. (|l996l ) treated for the first time 
the accretion rates of gas and solids in a self-consistent way. Their calculations 
show three distinct phases (fig. [1]): Phase I (t < 0.5 Myr) is characterized by a 
rapid build-up of the core. It ends when all the planetesimals in the core's initial 
feeding zone have been accreted. During phase II (0.5 < t < 7.5 Myr), the core 
is capable of extending its feeding zone by slowly accreting some surrounding 
gas. An increased core mass leads to a deeper potential well which leads to a 
contraction of the envelope which in turns leads to additional gas accretion from 
the surrounding disk and so on until the systems runs away and enters phase III 
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(i > 7.5 Myr). This phase starts at the critical mass which is reached when the 
mass of the core and the envelope become roughly equal. Runaway accretion 
occurs because in this regime the radiative losses from the envelope can no longer 
be compensated for by the accretional luminosity from the planetesimals. As a 
result, there is no equilibrium anymore and the envelope begins to contract on 
much shorter timescales. This contraction increases the gas accretion rate, which 
in turn increases the energy loses and the process runs away. The existence of 
such a critical mass is intrinsic to th e core-envelope m odel and does not depend 
upon the detail of the input physics (jStevensonlll982r). The critical core ma ss is 
usually of the order of 5 to 20 Earth masses ~(|Papaloizou Terauemlll999r ). In 
the runaway phase, the gas accretion rate is limited either by the planet itself 
(its ability to radiate away the gravitational energy) or by how much gas the 
disk can supply. 

The baseline core-accretion formation model has many appealing features, 
producing in a nebula with a surface density about four times the "minimum- 
mass" solar nebula a Jupiter like planet with an internal composition compatible 
to what is inferred from internal structure models (t hough the uncertainties in 
these models also allow a Jupiter without solid core (jSaumon fc Guillotl 12004^ . 
However, the timescale to form the planet (8 Myrs) has always been considered 
as uncomfortabl y long compared to observa tionally derived lifetimes of proto- 
planetary disks ( Haisch. Lada. Sz Ladal [200ll ) . Higher surface densities lead to 
significantly shorter formation timescales, which means that the baseline core 
accretion process is not intrinsically slow, but at the price tha t the resulting fi- 



nal content of heavy elements is very high (jPollack et al.| [l996V This "timescale 



problem" finally lead to the hypothesis that an other, faste r formation mecha- 
nism might be needed for giant gaseous planets (|Bosslll997l '). 



2.2. Improved and Extended Formation Models 



since iPollack et aD (|1996L the core accretion model has been significantly im- 
proved and extended, so that this "formation timescale proble m" is no longer a 
prob lem. Improved models use more accurate equations of state ( Saumon Sz Guillotl 
120041 ). more adequate solid accretion rat es dFortier et al. 1120071), take into ac- 
count local density patterns in the disk (jKlahr Bodenheimerl 120061 ) and the 
probable difference in opacity genera ted by grains in a p lanetary envelope as 
op posed to the interstellar medium ( Ikoma et al. I [2000). Indeed, as argued 
bv lPodolaki ([2003 ), grains entering the protoplanetary envelope may coagulate 
and settle out quickly into warmer regions where they are destroyed. The re- 
sulting lower opacity speeds up significantly the gas accretion timescale, while 
still fulfilling the obse rvational constraint on the abundance of heavy element 



(|Hubickvi et al. 1120051 ) 



Extended core accretion models have included new physical mechanisms: 
A concurrent calculation of the evolution of the protoplanetary disk, and most 
importantly, planetary migration. The discovery of numerous Hot Jupiters 
has forced upon us the necessity of planetary migration as in-situ formation 
of these objects is beyond the capability of any known formation theory. As 
planet formation, disk evolution and migration occur all on similar timescales, 
i t is necessary to treat these pr ocesses in a self consistent, coupled manner 
( Alibert. Mordasini. &; Benj[20 04). Extending the standard core accretion model 
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by t hese mechanisms n ot only leads to a natural explanation of the "Hot" plan- 
ets ( Alibert et al. 20061 ). but also solves the "timescale problem": For the same 
initial conditions that lead to runaway growth in the in situ case after 30 Myr, 
including concurrent disk evolution and migration leads to the formation of a 
Jupiter-like planet at 5.5 AU that has an internal composition compatible with 
inter nal structure models in just ~1 Myr starting with an embryo of 0.6 M® at 8 
AU (|Alibert et al.ll2005al ). The reason for this speed up is that o wing to migra- 



tion, t he planet's feeding zone is never as severely depleted as in lPollack et al 



(1996) and the lengthy phase II is skipped. Instead, the planet always migrates 
into regions of the disk where fresh planetesimals are available. 



3. Population Synthesis 



The improvements of the core accretion model in the past few years have made 
possible to account for the wealth of detailed in situ an d remote measurements 
that exist for Jupiter and Saturn ( Alibert et al. 2005bl ). showing that the core 



accretion paradigm has reached a significant level of maturity. In the same time, 
on the observational side immense progress has been made in both finding and 
characterizing new, extrasolar examples of the end products of the planetary 
formation process, but also the initial conditions for this process, i.e. the prop- 
erties of protoplanetary disks. These two developments have opened a new, fru it- 
ful p ossibility of bringing together observation and theory ( Ida Lin 2004al lbl. 
2005): The synthesis of populations of (extrasolar) planets from a core accretion 
model, which are then compared to the actually observed population. 

In this section we present some of the results of population synthesis cal- 
culations we have performed with our extended core accretion model. These 
simulations will be presented in a series of papers (Mordasini et al. in prep., 
Alibert et al. in prep.). Here, we can discuss only a very small fraction of our 
results, limiting ourselves to a host star mass of 1 M sun . 

3.1. Probability Distributions 

Population synthesis by means of a Monte Carlo approach rests on our ability 
to generate an appropriate set of initial conditions. To achieve this, the param- 
eters of the model specifying these initial conditions have to be drawn following 
suitable probability distributions. In our work, these probability distributions 
have been derived as much a possible from observations. 

The first Monte Carlo variable is the gas to dust ratio that we relate to the 
observed (stellar) meta llicity, using t he di stribution of [Fe/H] of the CORALIE 
planet search sample ( Santos et al. 2003). The second variable is the initial 
gas surface density, wh ich is constrained by the observed distribution of disk 
masses in p Ophiuchi ( Beckwith &: Sargentl 19961 ). The third variable is the 
photoevaporation rate which controls, together with the viscous evolution, the 
disk lifetime. We have adjusted the distribution of the photoevaporation rates 
so that the fraction of remainin g disk as a function of time is the same as 
observed by lHaisch et al. 1 ( 200ll ). The fourth Monte Carlo variable describes 
the initial position of the planetary seed we introduce in the protopl a netary disk. 
Since there are no observational constraints, we follow llda fc Lid (|2004al ) and 
distribute our seeds with an uniform distribution in log of the semi-major axis 
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in the regions of the disk where the local isolation mass exceeds the seeds mass 
(0.6 M®). Note that due to the relatively large seed mass and the other model 
restrictions, our synthetic planet population is incomplete in the mass domain 
of a few M©. In other words, our model is not a terrestrial planet formation 
model and should be completed appropriately before discussing this mass range. 

3.2. Detection Biases 

In order to compare to actual observations, once a synthetic population has been 
calculated we have to sort out the sub-population of planets which could have 
been detected by a given technique. For the radial velocity method, the detection 
probability depends to first order on the induced velocity amplitude and the 
instrumental accuracy. However, a large number of other quantities also affect 
the detection probability such as the magnitude of the star, its rotation rate, the 
actual measurement schedule, jitter, e tc. To com p ute th e actual observational 
detection bias, we use the method of iNaef et al.l ( 20051 ) which takes all these 
factors into account. 

3.3. Results 

Here we limit ourselves to the description of our nominal planet population 
model and restrict the host stars to solar-type stars. The most important pa- 
rameters characterizing the nominal case are a type I mig ration reduction facto r 
// = 0.01 and a disk viscosity parameter a = 0.01 (see lAlibert et al.l ( 2005al ) 



for an explanation of the quantities). To fix these parameters we compared the 
sub-population of detectable synthetic planets with known extrasolar planets or- 
biting single solar-type stars on low eccentricities orbits. We used Kolmogorov- 
Smirnov (KS) tests to assess the statistical signific ance of the synthetic mass- 
distance distribution (2D KS test. lPress et al.lll992l ). of the projected mass, the 



semimajor axis and the [Fe/H] distributions (ID KS tests), we also checked 
that we can reproduce the observed frequency of Hot Jupit ers as well as the 



stron g correlation of [Fe/H] and the detection probability ([Fischer &: Valenti 
2005). Fulfilling all these observational constraints simultaneously with a single 
synthetic population turned out to be a delicate task, as many combinations of 
parameters resulted in populations incompatible with the observational data. 

Mass-Distance Diagram The KS result for the two dimensional distribution in 
the mass-distance plane is the observational constraint that we weighted highest, 
as it is of similar importance t o planetary form ation as the Hertzsprung-Russell 
diagram for stellar evolution ( Ida Lmf|2004al ). In fig. [21 the left panel shows 



the a — M sin i of all synthetic planets. A close inspection of this figure shows 
that synthetic planets are not uniformly distributed but that some regions are 
more populated than others. The most obvious one is the vast population of 
failed cores, which are planets with masses of about 1 < M/M ffi < 10 seen at the 
lower mass boundary in the plot. At semimajor axes of about 0.8 to 4 AU and 
a mass of about 100 < M/M^ < 1000 a concentration of giant gaseous planets 
can be seen. This is the main clump as it corresponds to the location of giant 
planets for the most common initial condition. These planets typically s tarted 
at an initial semimajor axis just outside the iceline ( Ida &: Linl l2004al ) . For 



a < 5 AU, and 5 < M/M^ < 25, a horizontal bar is visible which corresponds to 
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Figure 2. Left: Minimum mass versus semimajor axis for all synthetic plan- 
ets of the nominal population. Right: The sub-population of synthetic planets 
detectable by a radial velocity survey with an instrumental accuracy of 10 m/s 
and duration of 10 years. Real extrasolar planets are indicated by large dots. 
Near the star subsequent evaporation of planets could be significant, making 
a direct comparison with the observations difficult. 

subcritical, core dominated, planets. These planets migrated inwards, mainly by 
type II migration, while only slowly accreting mass. Some of them grew massive 
enough to become super-critical thus leaving the horizontal bar towards the 
giant planet region. Others, never reached the critical mass and kept migrating 
inward building-up the population of Hot Neptunes. Above the horizontal bar, 
and to the left of the main clump, at masses between 25 to 200 M®, there is 
a region with some what less p l anets. This is the equivalent to the planetary 
desert described by llda Sz Lml (|2004al ). Compared to their calculations, the 
planetary desert in our calculations is certainly less pronounced. We note that 
this shallower planetary desert is not simply due to a sini effect. 

The right panel of fig. [2] shows the sub-population of synthetic planets de- 
tectable by a radial velocity survey with an instrumental accuracy of 10 m/s 
and a duration of 10 years. It is obvious that at that precision, only a small 
fraction (6.1%) of all synthetic planets can be detected. In the figure, the com- 
parison sample of actual extrasolar planets is also plotted. Comparing these two 
populations in a 2D KS test leads to a significance of 53% that both come from 
the same parent distribution. Comparing the Msini and the semimajor axis 
distribution separately in ID KS tests leads to a significance of 90% and 33%, 
respectively. The better KS result for M sin i than for a is attributed to the fact 
that in our models, migration is computed in a cruder approximation than mass 
accretion. 

Planetary IMF At a RV precision of 10 m/s, our knowledge of the planetary 
mass mass function was limited to the giant planets. At 1 m/s, the precision 
reached by HARPS, Neptune mass planets became detectable. In future, at 
a precision of 0.1 m/s, we will descend into the terrestrial mass regime. It is 
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Figure 3. Initial mass function of all planets of the synthetic population 
around G type stars. In the right panel, the population has been split in 
a low, medium and high metallicity bin. The region of a few Mq has been 
shaded as the model is incomplete there. 



therefore of interest to study the underlying, unbiased IMF (fig. [3]). From the left 
panel, it is seen that the IMF has a quite complex structure. Starting at the large 
mass end, we note that core accretion is able to form planets that can, at least 
if the presence of the core is not important, ignite deuterium burning. However, 
such planets are rare. At about 500 M ffi , the IMF has a local maximum, followed 
by a local minimum at ~ 40 M©. A small bump occurs in the Neptune mass 
domain due to the horizontal bar. At around 6 M®, the IMF finally starts to 
raise rapidly. As our model is incomplete for these very low masses, quantitative 
predictions should be regarded with caution here. Qualitatively, a strong raise 
of the IMF at such masses is however nevertheless expected, as it is simply a 
consequence of the fact that very often the conditions in the protoplanetary 
nebula are such that they don't allow the formation of a giant planet. 

In the right panel, the population was split in a low, medium and high 
metallicity bin. One can see that the IMF is clearly metallicity dependent, with 
metal rich systems producing more massive objects, and metal poor more small 
bodies. The distributions cross at around 6 M®. This metallicity dependent 
IMF explains why radial velocity technique based planet searches, which are 
biased towards large masses, have found planets preferentially orbiting metal 
rich stars. 

4. Conclusion 

The core accretion paradigm explains in an unified way the formation of giant 
and terrestrial planets, so that there is no need for a special mechanism for 
giant planets. Since the first core accretion models, significant improvement and 
extensions were made. Such improved and extended core accretion models can 
form giant planets well within observed disk lifetimes, so that there is no need for 
a faster formation mechanism. They have also reached a degree of maturity that 
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allows quantitative tests with observations, of both the giant planets of our own 
Solar System, and, by means of population synthesis, of the extrasolar planet 
population. In the latter case, the whole population of detected planets can be 
used to constrain the models, which excludes model fine tuning for a specific 
case, and fully exploits the observational investment. As shown by statistical 
tests, extended core accretion models can reproduce many observed properties 
and correlations in the extrasolar planet population in a quantitative significant 
way with one synthetic population at one time. This means that accretion 
models can now be used to predict future observations, so that theory can feed 
back on the design of future instruments 

Acknowledgments. We thank Isabelle Baraffe, Gilles Chabrier and Stephane 
Udry for useful discussions. This work was supported in part by the Swiss Na- 
tional Science Foundation. Computations were made on the ISIS and ISISII 
clusters at the University of Bern. 

References 

Alibert, Y., Mordasini, C, & Benz, W. 2005, A&A, 417, L25 

Alibert, Y., Mordasini, C, Benz, W., & Winisdoerffer, C. 2005, A&A, 434, 343 

Alibert, Y., Mousis, O., Mordasini, C, & Benz, W. 2005, ApJ, 626, L57 

Alibert, Y., Baraffe, I., Benz, W., Chabrier, G., Mordasini, C, Lovis, C, Mayor, M., 

Pepc, F., Bouchy, F, Queloz, D., & Udry, S. 2006, A&A, 455, L25 
Beckwith, S. V. W, & Sargent, A. I. 1996, Nat, 383, 139 
Bodenheimer, P., & Pollack, J. P. 1986, Icarus, 67, 391 
Boss, A. P. 1997, Science, 276, 1836 
Fischer, D. A., & Valenti, J. 2005, ApJ, 622, 1102 
Fortier, A., Benvenuto, O. G., & Brunini, A. 2007, A&A, 473, 311 
Haisch, K. E., Jr., Lada, E. A., & Lada, C. J. 2001, ApJ, 553, L153 
Hubickyj, O., Bodenheimer, P., & Lissauer, J. J. 2005, Icarus, 179, 415 
Ida, S., & Lin, D. N. C. 2004, ApJ, 604, 388 
Ida, S., & Lin, D. N. C. 2004, ApJ, 616, 567 
Ida, S., & Lin, D. N. C. 2005, ApJ, 626, 1045 
Ikoma, M., Nakazawa, K., & Emori, H. 2000, ApJ, 537, 1013 

Johansen, A., Oishi, J. S., Low, M.-M. M., Klahr, H., Henning, T., & Youdin, A. 2007, 

Nat, 448, 1022 
Klahr, H., & Bodenheimer, P. 2006, ApJ, 639, 432 

Mizuno, H., Nakazawa, K., & Hayashi, C. 1978, Prog. Th. Phys., 60, 3, 699 

Naef, D., Mayor, M., Beuzit, J.-L., Perrier, C., Queloz, D., Sivan, J. -P., & Udry, S. 

2005, 13th Cambridge Workshop on Cool Stars, Stellar Systems and the Sun, 

560, 833 

Papaloizou, J. C. B., & Terquem, C. 1999, ApJ, 521, 823 
Perri, F., & Cameron, A. G. W. 1974, Icarus, 22, 416 
Podolak, M. 2003, Icarus, 165, 428 

Pollack, J. B., Hubickyj, O., Bodenheimer, P., Lissauer, J. J., Podolak, M., & Green- 

zweig, Y. 1996, Icarus, 124, 62 
Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, 2nd cd. 

(Cambridge: University Press) 
Safronov, V. S. 1969 in Evolution of the Protoplanetary cloud and Formation of the 

Earth and Planets (Moscow: Nauka) 
Santos, N. C, Israelian, G., Mayor, M., Rebolo, R., & Udry, S. 2003, A&A, 398, 363 
Saumon, D., & Guillot, T. 2004, ApJ, 609, 1170 
Stevenson, D. J. 1982, Planetary and Space Science, 30, 755 
Wetherill, G. W., & Stewart, G. R. 1989, Icarus, 77, 330 



